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Abstract 

J> . The shape of a liquid surface bounded by an acute or obtuse planar angular sector 

is considered by using classical analysis methods. For acute angular sectors the two 
principal curvatures are of the order of the (fixed) mean curvature. But for obtuse 



■ sectors, the principal curvatures both diverge as the vertex is approached. The power- 

law divergence becomes stronger with increasing opening angle. Possible implications 
' of this contrasting behavior are suggested. 



PACS: 68.03.Cd — Surface tension and related phenomena 
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C : 1 Introduction 

o 
o 

The shape of the surface of a liquid or soap film constrained at its boundaries is a classic 
subject of mathematical physics fl], ||, [|. These studies demonstrate the power of producing 
[ very smooth surfaces of controlled curvature by choosing the shape of the boundary. The 

chief emphasis of these prior studies is to determine the global shape of the surface bounded 
by a given smooth surface. Here we emphasize the complementary question of the local 
surface shape in response to a singular boundary shape: namely, a line with a sharp bend 
enclosing a droplet spanning a plane sector of angle a. The role of singularities in governing 
the shape and the motion of fluids has aroused great current interest. Such singularities occur 
when a fluid droplet breaks apart [|, |, |, 0], when it merges another droplet || §, f7|, when 
it moves across a surface [|10l O, [L2|, O, ITJL ITS], [T7[, or when it moves through another 



fluid m 0, gg, |. 

Surprisingly, a qualitative change in the surface shape occurs as the opening angle of the 
boundary a increases past a right angle, as we show below. The curvatures for acute angles 
remain finite for the region near the vertex. But for obtuse angles, the curvatures diverge as 
the vertex is approached, with a power law that varies continuously with the angle. Similar 
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characteristic angles (not necessarily equal to 90 degrees) are encountered in the problem of 
capillary rise in a vertical wedge-shaped container. Here the meniscus height is bounded if 
the opening angle of the wedge is larger than some critical value and diverges as r — > if 
the opening angle is less than that critical value The contrasting behavior of acute 

and obtuse angles has also been noted for other phenomena involving Laplacian fields. In 
hydrodynamic flow, the velocity field near a wedge changes qualitatively as the angle increases 
through a right angle ||22|| . In diffusion (and analogous random- walk polymers) emanating 
into a wedge-shaped region there is a similar qualitative change of behavior [^3], ^ ^3] . 

Our motivation for focusing on droplets over an angular sector arises from observations of 
irregular droplets seen in everyday life. These often have sector-shaped regions arising from 
the vagaries of deposition and substrate shape. We have noticed that evaporation in these 
regions leads to distinctive drying patterns of solids dissolved in the liquid. To understand 
the nature of these drying patterns requires knowledge of the surface shape. For the circular 
drops the problem of the surface shape assumes very simple solution (spherical cap), allowing 
one to proceed with the issue of evaporation profiles up to the level of successful comparison 
of the theoretical results with the experimental data. These so-called "coffee-drop deposits" 
have aroused recent interest [P| ^7], E3. 

Specifically, we consider a droplet on the horizontal surface bounded by an angle a in 
the plane of the substrate (Fig. [[]). We assume that the droplet is sufficiently small so 
that the surface tension is dominant, and the gravitational effects can be safely neglected 
(significance of gravity increases with the size of the drop). At the same time, we do not 
assume that the contact angle between the "liquid-gas" surface and the plane is constant 
along the boundary line on the substrate. To achieve an angular boundary, the substrate 
must have scratches, grooves or other inhomogeneities (sufficiently small comparing to the 
dimensions of the droplet), which pin the contact line. A strongly pinned contact line can 
sustain a wide range of contact angles; the angle is not fixed by the interfacial tensions as it 
is on a uniform surface (Fig. Q). 

In the following section we first give a simple account of the shape that assumes that the 
liquid surface is nearly horizontal, and then we make a more systematic asymptotic analysis 
of the region where the distance r from the vertex is much smaller than the fixed inverse mean 
curvature R of the droplet, not making any a priori assumptions about the horizontalness 
of the surface. In the discussion section that follows, we calculate the curvatures of the 
obtained surface shape and describe some possible implications. In particular we discuss how 
the refraction of the light in the drop shows contrasting properties in acute versus obtuse 
angular regions. 

2 Calculation 

Boundary problem. Our purpose is to calculate the shape of the surface of the drop 
z(r,<p). Use of the cylindrical coordinates looks most natural in this problem so that the 
angle occupied by the liquid on the substrate is < r < oo and —a/2 < < a/2 and hence 
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(b) 

Figure 1: (a) A water droplet with a sector-shaped boundary on the plane substrate (side 
view), (b) The same droplet pictured from another point (top view; the experimental setup 
is sketched in Fig. 0). Black lines are the grooves on the substrate necessary to "pin" the 
contact line. (Photos by Itai Cohen.) 
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Substrate 



Figure 2: Illustration of the possibility of a wide range of contact angles in the presence of a 
groove or another inhomogeneity. 



the boundary conditions are 

z(0,<f)) = z{r,-a/2) = z{r,a/2) = (1) 
We start from the Laplace equation for the surface tension: 

2H=-^ (2) 
a 

where H is the mean curvature of the surface H = (ci + Qj)/2 and ci and C2 are the two 
principal curvatures. Here cr is the surface tension and Ap is the pressure difference between 
liquid and gas (Ap > 0). Since we neglect the effects of gravity, the pressure within the 
droplet is constant, and so is Ap. Thus Ap/a is just a constant parameter of the dimensions 
of inverse length. Since this is the only dimensional scale in the problem, introduction of the 
dimensionless variables r/R — > r and z/R —> z (where R = a/ Ap) makes the mathematical 
formulation parameterless: 

2H = -1 (3) 

Thereby we agree to measure all quantities of the dimension of length in units of R. Having 
found z(r, 0), one can restore the desired surface shape in ordinary units by simple substitu- 
tion z(r, 0) — > Rz (r/R, 0). 

Given a surface z(r, 0), its mean curvature can be computed in terms of the coefficients 
of the first and the second fundamental quadratic forms of the surface: 

TT I EN -2FM + GL 

H =2 EG-F* (4) 

Here E, F, G are the coefficients of the first fundamental quadratic form of the surface and L, 
M, N are those of the second fundamental quadratic form (see |29j or [|1|] for a concise review 
of these results of the differential geometry). For the surface x(r, 0) = r cos 0, y(r, 0) = r sin 
and arbitrary z(r, 0) these coefficients are: 

E = 1 + z 2 r F = z rZ(j) G = r 2 + zl (5) 
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and 



^ _ z _[z_ j^j- _ Zr <f> z 4>l r _ Z H> + rz r /„\ 

y/l + z 2 + (z /r) 2 + z r 2 + (^/r) 2 + z? + (z /r) 2 

Combining equations ©-© gives a second-order non-linear differential equation for the 
function z(r, d>): 



+ rz r + r 2 z rr + {z H zl + rz z r + z rr z\ - 2z r4> z r z 4 , + 2z r z\jr^ 



+ 



+ r 2 



+ {4 + (z,/rf}} 3/2 = (7) 
Note that this equation could also have been obtained if we had tackled the problem by 



minimizing the surface area A = JJ Jl + z 2 + (z^/r) 2 rdrdcj) (and hence the surface energy 
a A) while keeping the volume of the liquid beneath the surface V = JJ z rdrdcj) fixed[} This is 
equivalent to the minimization of the functional A — XV with respect to arbitrary variations 
of z(r,4>) that leave the boundary fixed, with A being a Lagrange multiplier. The proper 
choice of this parameter is A = Ap/a in ordinary units or A = 1 in dimensionless ones, which 
arises from the expression for the total energy E = a A — VAp. The Euler-Lagrange equation 
for the functional A — XV reads exactly as eq. (0). 

Thus, the boundary problem for z(r, <fi) consists of the equation (|7]) and boundary condi- 
tions (0). Note that we do not specify the boundary conditions at the opposite side of the 
drop (the furthest from the vertex), and this will lead to a set of undetermined coefficients 
in the solution to our problem. However, our purpose is to infer the universal features of 
the solution near the vertex, determined solely by the opening angle of the sector of interest 
and independent of the shape of the boundary outside of that sector. As we show below, 
knowledge of this subset of boundary conditions imposes sufficiently strict limitations on 
possible solutions, so that many important properties of the surface shape can be determined 
on the basis of only these local conditions. Had we specified all the boundary conditions, we 
would have obtained the exact parameter-free solution, dependent on the global shape of the 
boundary. 



Horizontal solution. There is no generic method for solving second-order non-linear 
differential equations of the kind of eq. ([?]), so we seek an approximate solution. First of all, 
we notice that if all partial derivatives of z are small (\z r \ 1 and \z$/r\ 1), i.e. if the 
surface is nearly horizontal, then the curly brackets in each pair of the square brackets can be 
neglected with respect to the rest of the terms. This horizontal approximation is not entirely 
obvious, and it will be justified in the next subsection. Thus, omitting the curly brackets in 
eq. (0), an easy-to-solve Poisson equation is recovered: 

V 2 z = -1 (8) 
lr The integrations are over the angular region occupied by the drop. 
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The general solution to the boundary problem (|S[), ([I]) can be written as a sum of three terms: 



z 



■— + z PN + Zgh 



where (— r 2 /4) is a solution to the non-homogeneous equation: 

/ ^.2\ 
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-1 



(9) 



(10) 



zpn is a particular solution to the homogeneous equation with non-homogeneous boundary 
conditions: 



V 2 z PN = and z PN (r, -a/2) = z PN (r, a/2) = — 



and Zgh is the general solution to the fully homogeneous boundary problem: 

V 2 z GH = and z GH (r, -a/2) = z GH (r,a/2) =0 (12) 
A particular solution to problem (ITT]) is 



r 2 cos 2<j> 
4 cos a 



Re( 



4 cos a 



if OL±\ 



z PN {r, 



^ In r cos 2(f) + ^0 sin 20 = Re (- ^) 



(13) 



a 



where a complex variable £ = re 1 ^ has been introduced on the plane of the substrate. Since 
this expression is a real part of an analytical function of £ (for each fixed a), it is a harmonic 
function on (r, 0)-plane by the Cauchy-Riemann conditions, and hence it is a solution to 
V 2 zpn = 0. Boundary conditions can be verified by direct substitution = ±a/2. 

The general solution to the homogeneous problem flT2| ) must satisfy the symmetry of the 
problem (i.e. must be even in <ft) and can be found by standard methods of mathematical 
physics: 



zgh 



n=0 



(2n+l)ir/a 



COS 



. 71(f) 

{2n + l — 

a 




(2n+l)ir/a 



(14) 



The constants C n cannot be determined without imposing further conditions on the solution 
(for instance, obviously C n may depend on a). Had we specified the boundary conditions 
along some curve r o (0) that represents the rest of the boundary of the drop, all the C n would 
be fixed. Since the number of coefficients C n is infinite, any reasonable boundary condition 
at ro(</>) can be satisfied. On the other hand, as it is apparent from our construction, those 
missing boundary conditions would not influence any other terms in the solution, which are 
universal and do not depend on the rest of the drop. 

Thus, the general solution (even in and going to as r — > in cylindrical coordinates) 
to the boundary problem (H), ([!]) is 



z(r, 0) 



+ r 4 CJ ^ + E C n r^ 2n+l >l a cos (2n + IV 

4 4 cose? " V ' 

n=0 



if a±\ 



^ + H0sin20 



lnr • cos 20) + £ C n r^ 2n+1 ^ 2 cos [(2n + 1)20] 

n=0 



if 



a 



(15) 
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Note the dominant terms in the limit r <C 1 for different values of a: for acute angles 
(a < it/2) the r 2 -term dominates, for obtuse angles (a > tt/2) the W a -term does, and for 
the right angle (a = vr/2) the (r 2 lnr)-term does. At a = tt/2 both the r 2 and the W a terms 
scale with r as r 2 (i.e. they "switch" here in the sense of power dominance), and at exactly 
this value a logarithmic correction to r 2 appears, as it typically happens for a power series 
solution near a crossover of two powers. 

It may seem from the structure of the expression (|l"5|) that this solution is a discontinuous 
function of a at a = 7r/2. However, this is not true. The key observation is that the 
coefficients C n can be different for different values of a. In particular, Co in the upper line 
of eq. (|lf|) is not the same as the one in the lower line. Let us keep the notation Co for the 
coefficient Co in the right-angle expression (lower line) and introduce a new notation C for 
that coefficient in the expression for angles different from tt/2 (upper line). Consider some 
angle a in the vicinity of tt/2, i.e. let a = tt/2 + e, where \e\ <C 1, and expand the second 
and the CW" cos(TT(fr/a) terms in the result for a ^ tt/2 in small parameter e: 

r 2 cos 20 7r/o! t ± i \ r2 ( cos 20 



4 cos a 



+ CW" cos(7r0/a) = '- (l + 0(e 2 )) + 



+ Cr 2 (l - ^-e In r + 0{e 2 )^j (cos 20 + ^e0 sin 20 + 0(e 2 )^j (16) 
Now, since C can depend on a (and hence e), we set 

C = ±- + C + O(e) (17) 

Then the two diverging terms of the order of 1/e cancel, and we recover (up to the leading 
order in e) the second and the Cor n ^ a cos(7T0/o;) terms in the result for a = n/2: 

r cos20 + (7r 7r/a/ , / n = _!L i n r cos 2 + ^0 s in 20 + C r 2 cos 20 + 0(e) (18) 

4 COS a 7T 7T 

Since the other terms in the result ([15]) are identical in the upper and the lower lines, we 
have shown thereby that our solution is indeed continuous in a at fixed r for a = tt/2 (or 
e = 0). 

Thus, solution (|T3|) behaves reasonably well in the full range of values of angle a from to 
tt. It cannot be used for a > tt since in that range it violates the horizontalness requirement 
employed in its derivation. 



Asymptotic analysis. The results above required the assumption that the drop is 
nearly horizontal. This assumption has not been justified yet, and now we justify it via a 
more systematic treatment. Since we are interested in the behavior of the surface near the 
vertex of the angle, we introduce a new small parameter for the problem 

r ^ R = — — (ordinary units) (19) 
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or 

r < 1 (dimensionless units) (20) 

For small r we may write z(r, <fi) as a standard series expansion: 

z (r,<p) = r ^(0)+ r ^(0) + --- (21) 

where < v <//<.. .. Note that we do not restrict our attention to the horizontal case only, 
i.e. we do not require 1 < v. Values of v between and 1 leading to non-horizontal surfaces 
will be eliminated automatically by application of boundary conditions to the solutions of 
ec L- ©j thus justifying the horizontal assumption. Here we find only the main asymptotic 
(z/-term) and the first order correction (/x-term), but the method allows one to proceed up 
to an arbitrary order. Details of the calculation are considered in the Appendix; results 
are presented below. We also treat the case of the right angle separately since we expect 



logarithmic corrections to the main power of r and failure of the assumption (El]) 



Leading asymptotic. Substitution of z(r,<f>) = r u <& v ((f)) into eq. (|7|) and retention of 
only the dominant terms for r C 1 lead to different equations for different possible values of 
power v. Solution of those equations and application of symmetry arguments and boundary 
conditions eliminate someQ values of v, leaving at the end only two possibilities iy = 2 and 
v = ii /a). For these two values the terms retained are a subset of those constituting eq. (§), 
yielding the following main order result: 

f if 0<«<f (u = 2) 

z{r,<f>) = l (22) 
[ CW Q cos ^ if f < a < 7T (y = it /a) 

This agrees with the leading behavior of the horizontal solution fll5|) as r —>■ 0. Thus, our 
surface is indeed nearly horizontal (since v > 1), and the horizontal approach indeed produced 
a sensible result. 

The constant C is again restricted by neither the equation nor the side boundary condi- 
tions, but it would get fixed once the boundary conditions at the furthest side are taken into 
account. It is a direct equivalent to the constant C in the horizontal solution fll5|) obtained 
by an independent treatment (recall that we relabeled the Co in the upper line of eq. ( |T5| ) 
into C). 

Dependence v(a) is shown on Fig. |3|: for acute angles v = 2 and for obtuse ones v = tt/cx. 

First order correction. The final first-order result for the surface shape in the limit 
r <^ 1 is: 

I. If < a < 7r/4, then v = 2, // = 4 and 

z (r, 4>) = -t 1 + 

4 y cos a J 

2 The leading-order angular solution for v < 1 is = Ccos^^. This cannot vanish as required for 

4> —> ±a/2 (see Appendix for details). 
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0- 

Jt/4 it/2 3n/4 

Angle a 

Figure 3: Dependences z/(a) and //(a). 



cos 40 



2 tan 2 a + 12 



+ r 4 4tan 2 a-3 r + 

192 \L J cos 2a 

2. If 7r/4 < a < 7r/2, then v = 2, = 7r/o; and 



r 2 f CQS2 ^ _ 1 



cos a 



cos 20 
cos a 



6 tan a + 9 



C W a cos^ + - 



a 



(23) 



(24) 



3. If 7r/2 < a < 37r/4, then v = n/a, fi = 2 and 



*(r, 0) 



<>*/« cos ^ + -r 2 f 
a 4 \ 



cos 20 



1 + 



cos a 



4. If 37r/4 < a < it, then v — ir/a, fi = 37i/a — 2 and 

r 3 

+ ^ 3 — 77- 



z r, 



Cr* la cos — + 



4a 2 (27r - a 



. r (*r/a)-2 00S !^ + 



(25) 



(26) 



The structure of the solution (p3[)-(p6|) becomes clear if one plots functions via) and //(a) 
together (Fig. [3|). Four powers of r appear in these formulas: r 2 , r 4 , W Q and r^ 3n ^ a ^~ 2 . For 
any given a our procedure selects the two lowest powers in this set of four. Different powers 
get selected for different a; this leads to the four cases appearing in (|23| ) - (|26|) . In a full 
expansion, we expect all four powers to be present for all angles^. 

Since a plane z = Dx = Dr cos is the exact solution to the equation (|7|) when the 
boundary of the surface is a straight line a — 7T, it is rewarding to observe that both terms 
in the expression (p6~[) reduce to this functional form with D = C + (C 3 /4) + • • • as a — > tt. 

3 Note that all terms satisfy boundary conditions independently. 
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Note that although the leading asymptotics in r is the same in the results flI5p and fl2"3p- 
(p6|), the sub- leading terms are different. This is due to the fact that the two results are based 
on different approximations: the former assumes \z r \ <C 1 and \z$jr\ <C 1 while the latter 
assumes r C 1. Since, as shown in the previous subsection, the horizontal approximation 
follows from the close-to-the- vertex one (r <C 1), the asymptotic treatment of this subsection 
describes the surface shape more accurately than result (0), picking up lower powers of r 
for the first sub-leading terms. 

Right-angle sector: first two terms in the expansion. We already know that a 
pure power series does not work in the case a = tt/2 and that the leading power of r should be 
close to 2, at most logarithmically close. So, we introduce a new ansatz instead of series fl2T|): 

z(r,<f>) = (-r 2 lnr)$ 1 (0)+r 2 $ 2 (0) + --- (27) 

Subsequently, by repeating the steps of the Appendix, we find $i(0) = v4cos20, then obtain 
the second term ($ 2 (0)) and fix A = 1/tt by boundary conditions. As a result we recover 
the lower line of the expression (|T5|): 

z(r, 0) = --r 2 In r cos 20 + r 2 f-0 sin 20 - - + C cos 2<p) H (28) 

7T \7T 4 / 

where Cq is again a constant equivalent to the Cq in fll5l). The relation between the constant 
C in this right-angle expression and the constant C in the result (|2"B|)-([2"ED for a ^ tt/2 is 
exactly the same as in eq. (|17D due to the continuity of the full solution for all values of a. 

Thus, an asymptotic expansion of the dimensionless function z of the dimensionless co- 
ordinates r and in the limit r -C 1 has been found for all values of a in the range from 
to 7i. A natural extension of this study is to consider angular regions that have a > tt. We 
guess that for such cases z(r, 0) grows as a power of r less than unity, so that the slope of 
the surface diverges at the vertex. But we have not succeeded to verify this behavior with 
our methods. 

Typical behavior of the universal function z(r,<p) (eqs. (P^[)-(P^) and fl2"8])) for angles 
a = 37r/8, 7i / 2 and 5tt/8 is shown on Fig. ^. To facilitate the comparison, we also plot 
the bisector cross-sections of the same three surface profiles in one frame (Fig. ^). Constant 
Cq is taken to be 1 in all cases, while the value of constant C is chosen according to the 
prescription (|17|). In spite of the existence of drastic mathematical difference between the 
three regimes, this difference is not apparent from looking at the profile of the surface itself 
(as on Fig. |j), and that is probably why it went unnoticed so far. We will further emphasize 
this qualitative difference in the following section. 

3 Discussion 

Curvature. Let us better understand the main order result by looking at the curvature 
of the surfaces described above. The principal curvatures along the bisector of the angular 



10 



o.: 

(a) 




Figure 4: Surface shape z(r, <p) (first two terms of the expansion in r with (7 = 1/ (4a— 2tc)+Cq 
and C = 1) for (a) a = 3tt/8 (eq. (b) a = vr/2 (eq. @), and (c) a = 5vr/8 (eq. (^)). 

The mathematical differences in the shapes are not apparent in this view. 



11 





a = 371/8 




0.008 - 


a = nil 

a = 571/8 




0.006- 




0.004- 




0.002- 


/ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 


0.000- 





0.00 0.02 0.04 0.06 0.08 0.10 

Distance from the vertex r 

Figure 5: Surface profiles at bisector (0 = 0) for three values of angle a. 



region are in the f and directions, and for the points on the bisector the radial (c rr ) and 
the azimuthal (c^) curvatures are simplyf] 



L 
E 



and 



N 
G 



(29) 



with E, G, L and N defined in eqs. (^)-@. For small r (and therefore small \z r \ and |z</j/ r l) 
these expressions simplify even further to 



and 



(30) 



so that the principal curvatures sum to V 2 z. Thus, for the surface fl22l) and (^8|) (up to the 
leading order in r 1 only) these principal curvatures are: 



■ 1 + ^— 

2 1 2 cos a 



1 



1 



c rr = — - In r 



£, n(ir-a) r ( n / a )-2 



2 2 cos a 

- lnr 

7T 

^ 7r(7r-a) ^( 7r / a )-2 



if < a < tt/2 

if a = ix j2 (31) 

if tt/2 < a < 7r 



Typical behavior of radial curvature c rr on the bisector as a function of r is shown on Fig. ^| 
for the same three values of a as on Figs. [| and |5|. Now the dramatic difference between the 



4 For an arbitrary point of the surface the principal curvatures are determined by the roots c\ and C2 of 
the quadratic equation 

T-rK M — of? 

= 



L-cE M - cF 
M -cF N ~cG 



At the points where the principal directions are along the coordinate ones (e.g. on bisector) coefficients F 
and M vanish and the principal curvatures are given by simple relations (p9|) (see [fil p9| for details). 
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Figure 6: Radial curvature on the bisector for three values of angle a (main asymptotic only). 



two regimes of a becomes apparent. For acute angles the curvatures remain finite as r — > 
while for obtuse ones the curvatures diverge as a negative power of r (changing from to — 1 
as angle a passes from 7r/2 to 7r). The limiting case of the right angle has an intermediate 
logarithmic divergence. 

Note that the finite values of curvature (for a < 7r/2) sum to —1 in full accord with 
equation (|). However, the divergent values (for a > 7r/2) have the opposite signs and thus 
sum to 0. This is a result of the neglect of the corrections to the main asymptotic in r. 
Had we kept the corrections to the divergent curvatures, they would sum to —1. Thus, for 
instance, the second term in the solution (|25| ) is the kind of the correction which provides 
summation to —1 when the main (divergent) terms sum to 0. 

The origin of the arbitrary constant C in the solution (|22|) is now seen to be related to the 
divergence of principal curvatures for obtuse angles. Indeed, for acute angles the curvatures 
are finite and c rr + c<^ = — 1 while for the obtuse ones the curvatures are divergent and 
c rr + = (up to the main order in r). Therefore, multiplication of the solution z(r,<f>) 
by an arbitrary constant is not allowed in the former case while is perfectly legitimate in the 
latter one (since both c rr and get multiplied by the same constant). Thus, a possibility 
for an arbitrary multiplicative constant in the solution for obtuse angles comes from the 
divergence of the curvature, which in its turn reflects a different r-dependence of z(r, <fi) for 
a > 7r/2. 

On the other hand, the presence of an undetermined coefficient means that the shape 
is influenced by the boundary conditions at the side of the drop furthest from the vertex. 
Hence, different scaling with r for different angles results in different dependence on boundary 
conditions outside of the sector of interest: for acute angles the dimensionless shape of the 
surface at the tip of the sector does not depend on these boundary conditions while for 
obtuse angles it does. Of course, this argument was based on the main-order solution (|22|), 
but general dependence of the surface shape on the boundary conditions at the opposite side 
of the drop via constant C is already apparent from the first-order solution (|23|)- (|26|) : the 
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Figure 7: Experimental setup for the refraction demonstration. 



larger the opening angle, the stronger dependence of z(r, <p) on these yet-unspecified boundary 
conditions (the number of terms containing C increases as a increases). This seems quite 
reasonable, as intuitively surface shape near the vertex must cease depending on the rest of 
the drop as a — > 0, and it must be fully specified by the rest of the drop when there is no 
vertex at all (i.e. when a = ir). 

Experimental realizations. The contrast between the drops over acute and obtuse 
angular regions may be seen in the way they refract light. To illustrate, we picture a pair 
of sector-shaped droplets on a transparent substrate at a distance s above an object plane 
(Fig. |7|). One angle is acute, and the other one is obtuse. The object plane consists of a set of 
closely spaced parallel lines perpendicular to the bisectors of each angular region, so that the 
spacing between the lines is in the f direction along the bisector. Observation of the object 
plane through the droplets allows one to make qualitative judgment about the behavior of 
the curvature near the vertices of each sector. The result of such a simple demonstration 
with spacing between the parallel lines of approximately 1.6 mm, acute angle of about 51°, 
and obtuse angle of about 124° is shown in Figure |j. As it is apparent from this image, the 
spacing between the lines seen through the drop over the obtuse angular region decreases as 
they approach the vertex, while the spacing between the lines seen through the drop over 
the acute sector remains unchanged. Note that only a few millimeters near the vertex should 
be taken into account while viewing this figure since the inverse mean curvature R for water 
drops not distorted by gravity is of this order of magnitude. At higher distances gravitational 
effects cannot be neglected while calculating the surface shape. 

This observation agrees favorably to the result of our calculation, as it can be seen from an 
argument based on geometrical optics. Indeed, for the dimensions of the optical image along 
the bisector only radial curvature is important, and one can write the following approximate 
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Figure 8: A set of parallel lines as seen through the drops over acute (left) and obtuse (right) 
angular regions. The spacing between the undistorted lines is about 1.6 mm. The opening 
angles of the acute and obtuse sectors are approximately 51° and 124° respectively. The 
regions of interest are only a few lines (a few millimeters) around the vertex of each drop. 
The grooves necessary to pin the contact line can also be seen on this image. (Photo by Itai 
Cohen.) 

expression for the linear magnification by the drop in the f direction: 



1 + s(n — l)c rr 

where n is the index of refraction of the liquid the droplets are made olQ. According to our 
result, for obtuse angles the curvature diverges and the magnification should go to zero as 
the vertex is approached; for acute angles both quantities remain finite. Thus, qualitative 
validity of our result is confirmed by the simple demonstration described above. Similar 
behavior should hold for the light reflected off the surface of the droplet because curvature 
is equally important for both such phenomena. 

Another possible system to test the predictions of our study is a pressurized soap film. As 
shown in Figure |9|, a soap film on a wedge-shaped frame with an applied constant pressure 
difference across it will have constant mean curvature. Thus, it is described by our formalism. 
The only difference between such a film and a liquid drop is that the film has two surfaces, 
and therefore the applied Ap must be two times as much as Ap in equation @. Once such 
a surface is produced, it can be made as big as necessary for experimental convenience, since 
gravitational effects are virtually absent for this realization. 

As these examples illustrate, the change in behavior on going from acute to obtuse planar 
angular regions can show up in concrete ways. We suspect that further differences will 
emerge as capillary flow and evaporation properties of these sector-shaped liquid interfaces 
are explored. 

5 This expression assumes the horizontalness of the drop. 



15 



Soap film 
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Appendix 

Here we present some details on how we obtain expressions (22) and (^)-(|26|). One starts 
from the substitution of z(r, <p) = r iy $ i/ (0) (the first term of the expansion fl2~ID) into eq. (|7|) 
and obtains the following equation for Q u (<p): 



(P[Q V ) + r 2 ^ 2 Q[^]) + r 2 (l + r 2 - 2 R{^ u }f 2 = 



where 



P[$ v ] = C + v 2 $ v 
Q[$ v ] = v 2 $ 2 M + - ^ 2 ) *v (Kf + v^l 
R[$ u ] = (Kf + v 2 <b 2 v 

Considering all possible values of u, leaving only main terms in r (the smallest powers of r) 
and solving for Q u (<f)) in each case, one arrives at the following set of solutions (only even 
terms are shown due to the symmetry of the problem): 



( Ccos" 



C COS0 



C cos vc 



C cos 20- \ 



no solution 



if < v < 1 

if u = 1 

if 1 < v < 2 

if v = 2 

if v > 2 



(C is independent of r and <ft everywhere but arbitrary otherwise.) Obviously, the first two 
options can not satisfy boundary conditions $^(— a/2) = $ y (a/2) = for angles a < n, and 
thus the "not horizontal" solutions with v < 1 are naturally eliminated. In cases 3 and 4 
boundary conditions yield v = n/a and C = 1/(4 cos a) respectively. Thus, the main order 
result in the limit r 1 is nothing but eq. (f22|): 



z(r, 0) 



cos 2(j> 



- 1 



CW« cos i± 



if < a < f 



if | < a < 7r 



(i/ = 2) 
[y = 7t/a) 
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Then we proceed in exactly the same fashion to determine \i and to find $ M (0) by employ- 
ing just calculated main order result. Substitution of the first two terms of the expansion (|2T|) 
into equation (0) yields the following equation for 

r u (P[$„] + r"-"^] + r 2 »- 2 Q[<$> v ] + r u+ ^ 2 V[^, $ M ] + 0(r 2 ^ 2 )) + 

+ r 2 (l + r 2u ~ 2 R[^ u } + r^~ 2 iy [$„, + 0{r 2 ^ 2 ]f /2 = 
where <& v {4>) is already known, P, Q, R are the same as in the equation for <& u ((j)), 

= + /i 2 $ M 

and expressions for and W are irrelevant to any final results. Similar to the previous 
case analysis, including thorough consideration of all possible cases for values of fi, neglect 
of the terms of the order higher than the first correction in r and application of symmetry 
arguments and boundary conditions to the solutions, leads to the first-order result (|23|)-(|26|). 

Obviously, the procedure of building the next term of power series for the solution of 
eq. (0) can be repeated up to an arbitrary order. 
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